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Gaussian states — or, more generally, Gaussian operators — play an important role in 
Quantum Optics and Quantum Information Science, both in discussions about concep- 
tual issues and in practical applications. We describe, in a tutorial manner, a systematic 
operator method for first characterizing such states and then investigating their prop- 
erties. The central numerical quantities are the covariance matrix that specifies the 
characteristic function of the state, and the closely related matrices associated with 
Wigner's and Glauber's phase space functions. For pedagogical reasons, we restrict the 
discussion to one-dimensional and two-dimensional Gaussian states, for which we provide 
illustrating and instructive examples. 
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1. Introduction 

Gaussian wave functions appeared already very early in the development of Quan- 
tum Mechanics, when Schrodigner explained the observed behavior of the linear 
harmonic oscillator in terms of the undulatory eigenfunctions of the correspond- 
ing differential equation.^ Heisenberg, in his fundamental work of 1927 devoted to 
the uncertainty relation,^ used a Gaussian wave function of the form (in historical 
notation) 

-iaiVf. _ 221// _ „') 
5(?7, q) prop e ^ ^' (1) 
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to show that the undulatory character of this wave function, combined with the 
Dirac-Jordan transformation theory, leads to the indeterminacy relations. Heisen- 
berg used this Gaussian wave function as a probability amplitude — perhaps the 
first application of this kind — to exhibit, in his own words, that The more pre- 
cisely the position is determined, the less precisely the momentum is known in this 
instant, and vice versa. 

Gaussian statistical properties arc fundamental to many applications in statisti- 
cal physics. It is familiar that the classical theory of a Gaussian noise associated with 
Brownian motion is fully characterized by the covariance function of the phase space 
variables q and p. The transition from classical phase space variables to canonical 
quantum position and momentum operators requires "quantization rules" for asso- 
ciating quantum operator functions F{q,p) with the classical numerical functions 



In the framework of quantum noise associated with a one-dimensional harmonic 
oscillator, Gardiner"^ was perhaps the first to discuss, in a textbook format, various 
phase space properties of the most general Gaussian operators. Such a systematic 
formalism is very useful for the discussion of Gaussian statistical properties of a 
single quantum system with one degree of freedom. 

The remarkable role that is played by one-dimensional Gaussian wave func- 
tions in realizing the ultimate limit of Heisenberg's uncertainty relation for a sin- 
gle position-momentum pair naturally invites generalizations to higher-dimensional 
systems. That, however, opened a path into the then-unexplored territory of quan- 
tum correlations and quantum entanglement and thus started a never ending story. 

In their reasoning concerning the alleged incompleteness of quantum 
mechanics,^ Einstein, Podolsky, and Rosen (EPR) used in 1935 the following wave 
function for a system composed of two particles (in historical notation): 



The fact that this function could be written as a sum of products of one factor each 
for the two particles. 



leads to the intriguing concept of quantum entanglement, a term coined by 
Schrodinger in Refs. 5 and 6.^ 



''In Heisenberg's more general context, this Gaussian wave function refers to an arbitrary one- 
dimensional system, be it a free particle, a harmonic oscillator, or something with other dynamical 
properties. For the sake of convenience, we shall invariably speak of harmonic oscillators, thinking 
in particular of those associated with modes of the quantized radiation field. 
''It was in Ref. 5 that the word entanglement, which refers to the non-separability of quantum 
states of composite systems, was given its familiar meaning within the context of quantum me- 
chanics. For a brief history about the use of the English word entanglement and the German word 
Verschrankung see Ref. 7. 



T{q,p). 
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The entangled EPR wave function (2) is a singular function of the distance 
x\ — X2 and can be visualized as an infinitely sharp two-dimensional Gaussian wave 
function of the entangled two-party system. Soon after the publication of the EPR 
paper, Bohr pointed out* that the two pairs of canonical variables qi,pi and 92, P2 
for the two particles of a composed system can be replaced by two new pairs of 
conjugate variables, 

Qi,2= -^{qiTq2), Pi,2 = -^{Pi T P2) , (4) 

each pair now referring to both particles. Bohr noted that, since they commute: 
21^2 = 'P2Q1, the obscrvablcs Qi and V2 can be assigned sharp values simultane- 
ously, so that the wave function (2) represents a joint eigenstate of these commuting 
variables. This property is at the heart of quantum entanglement, unintentionally 
brought to light by EPR. Bell inequalities of some kind are violated for the EPR 
wave function, as can be demonstrated by using its Wigncr representation.^ 

Schrodinger, Heisenberg, and the EPR trio were dealing with systems described 
by wave functions, i.e., with pure quantum states. For a quantum system described 
by a Gaussian mixed statistical operator p, the correspondence between classical 
Gaussian functions G{qi,Pi,q2,P2) in two dimensions and the statistical operators 
p = G{qi,pi,q2,P2) should employ a "quantization rule" that preserves all the 
properties of a density matrix. 

Recent applications of entangled two-mode squeezed states of light for quantum 
teleportation^° and other quantum information purposes^^ have generated a lot of 
interest in the entangled properties of general mixed Gaussian states in quantum 
optics. "'^^ It turns out that the concept of quantum entanglement, as defined by 
(3), has to be generalized when the system is not in a pure state. In the general 
case of a density operator, rather than a wave function, one uses the definition of 
quantum separability introduced by Werner: A general quantum density operator 
of a two-party system is separable if it is a convex sum of product states, 

p = ^2 ^fe Pi^^ P2^^ witli ^2 ~ Wk > 0, (5) 

k k 

where p^''^ and pj'^-' arc statistical operators of the two subsystems in question.*^ 

The separability properties of general mixed states of a harmonic oscillator in 
two dimensions can be studied with the help of two different methods. The first 
uses the covariance matrix and the Heisenberg uncertainty relations.^'* The second 
makes use of the criterion of positivity under partial transposition.^''' These two 
approaches lead to essentially the same conclusions, while employing very different 
techniques. 



'^Actually, one should only require that the given p can be approximated to any required accuracy 
by a sum of this kind with a finite number of terms, but we take the liberty to ignore mathematical 
details of such a more pedantic sort. 
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The objective of this tutorial is to review different operator and phase-space 
techniques that allows a systematic investigation of different properties of Gaus- 
sian operators of unit trace, referring, for instance, to a single mode or two modes 
of the radiation field. We discuss the properties of Gaussian operators with the 
aid of various techniques that are widely used in Quantum Optics. In particular, 
we provide a careful description of Gaussian operators using (i) the boson opera- 
tor algebra, and (ii) phase-space descriptions based on the Wigner representation 
and the Glauber P-representation. We establish relations between these equivalent 
though different-in-form versions of Gaussian statistical operators. In Sec. 2 we 
gather the necessary mathematical tools used in this tutorial at the example of the 
general Gaussian operator of a one-dimensional harmonic oscillator. Wc provide a 
general parametrization of such states, discuss the positivity criteria and the so- 
called P-representability of positive Gaussian operators. At the end of Sec. 2 we 
illustrate the formalism with the example of a squeezed single- mode state of light. 
Then, in Sec. 3, we turn to two-party Gaussian states. The central part of this tuto- 
rial is Sec. 3.4 that deals with the separability of two-dimensional Gaussian states. 
We discuss the positivity criteria of such states and their separability conditions, 
thereby illustrating various approaches. At the end of Sec. 3 we provide a number 
of simple examples that are useful for studying and illustrating various physical ef- 
fects in Quantum Optics^^ and in Quantum Information Theory. We confine our 
discussion to two-party systems, but the methods presented in this tutorial can be 
extended to many-party systems. In concluding remarks we summarize and suggest 
some supplementary reading on the subject. 

2. Gaussian States of a One-Dimensional Oscillator 
2.1. Parameterizations 

Any operator referring to a harmonic oscillator — position operator q, momentum 
operator p, both measured in natural units, so that qp — pq = i — is a function of 
the familiar ladder operators 

. q-ip q + ip 

We can specify such an operator G{a^ ,a) by its characteristic function C{z*,z), 

C{z\z) = Tr {e^"^ " ^*"G(at, a)} , (7) 
which is a numerical function of the complex phase space variables 

Here, q',p' are the cartesian coordinates of classical phase space as one knows them 
from Hamilton's approach to classical mechanics or the Liouville formulation of 
statistical mechanics. A more compact way of writing (7) is 

C(z) =TY{e-^^^"G(a)} (9) 
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where we introduce 2-component columns and rows in accordance with 



and meet the symplectic 2x2 matrix 

E = 



a = 



1 
-1 



a^ = [a\a] , 



(10) 



(11) 



Reciprocally, we get the operator from its characteristic function by a phase space 
integration, 



G(a)= J{dz)e-^^'C{z). 



(12) 



which is essentially the Liouville trace of statistical mechanics and complements the 
quantum- mechanical trace of (9). The explicit form of the volume element (dz) de- 
pends on the way phase space is parameterized. For the standard parameterization 

(8) one has 

dq' dp' 



(dz) = 



27r 



(13) 



or (dz) = ;i:dRe^ dlm^; = ^d^z in a popular notation. But sometimes the actual 
integrations are more easily performed with 

dq'dp' 



z = q 



z = —ip 



(dz) = 



27r 



(14) 



for instance. Note that we include the normalizing denominator of 27r in (dz) which, 
in a rough manner of speaking, indicates that the integral in (12) counts one quan- 
tum state per phase space volume of 27r=27rfi. 

We focus on Hermitian operators G(a) = [G(a)] ^ with characteristic functions 
that are of the Gaussian form 



C{z) 



-iztCz 



lmz*2 1)2*2 _ 1^*22 



with the 2x2 matrix C given by 



C = 



n + 
m* 



n + 



= ct. 



(15) 



(16) 



The hermiticity of G implies [C(z)]* = C(— z), and vice versa; the Gaussian char- 
acteristic function (15) has this symmetry property because the parameter n is 
real. 

For a given C(z), (15) does not specify the diagonal entries of C uniquely, only 
their sum is determined. The symmetry of the standard form (16) of C exploits this 
arbitrariness conveniently. A particularly useful way of looking at this proceeds 
from the response of a and z to transposition. 



= [a,a^] =a^T, z"^ = [0, z*] = z+T , 



(17) 



6 B.-G. Englert and K. Wodkiewicz 



with 



SO that 



(ztCz)^ = ztTCTTz. (19) 



This states that one can replace C by TC'^T in (15), or by any hnear average of 
both, without changing the characteristic function C(z). It is therefore natural to 
enforce the symmetry 

C = TC^T , (20) 

and this establishes (16). 

Since C(z) = 1 for z* = 0, z = in (15) wc take for granted that G has 
unit trace, Tr {G} = 1. This excludes Gaussian operators with infinite trace, but 
otherwise it is just a matter of convenient conventional normalization. 

The absence of linear terms in the exponent indicates another convention: We 
assume that Tr {aG} — 0, Tr|o'''G} = 0, which can always be arranged with a 
suitable unitary shift of and a. 

Upon expanding the characteristic function in powers of z* and z one easily 
identifies the physical significance of the numerical parameters n and m, 

n = Tr {a'faG} , 

m = -Tr {a'^G} , m* = -Tr {a^^c} , (21) 

or, more compactly, 

Tr{aa1^G} = ECE+ iE (22) 

and 

C= ETr{aatG}E-iE= iETr{[aa1' + T(aa1^)'^T]G}E. (23) 

These relations identify C as the covariance m,atrix of G. 

Anticipating that this will be of some relevance later, we note that a positive G 
can serve as a probability operator (alternatively called "state operator" or "density 
operator"). Then Heisenberg's uncertainty relation requires 

{a^a){aa)) > {a^^){a'') (24) 
so that n and m of a positive G must be such that 

n{n + l)>m*m or n > -y/ m*m + \ — \- (25) 
As a compact statement about the covariance matrix C, this appears as 

C + iE>0. (26) 
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Here we recall that one derivation of (24) simply exploits the positivity of operators 
of the form X'^X. Choose 



X = a^Ex with 

so that 



(27) 



(x+Eaa^Ex) = xtE(aa^)Ex > (28) 

has to hold for any x. With (22), the recognition that (for G > 0) 

E(aa^)E = ETr{aa+G}E = C + iE (29) 

then establishes (26) as a necessary property of any positive G. 

In (12), G is expanded in the Wcyl basis'^ that consists of the unitary operators 
e^ = e~^ Equivalently, we can use the Hermitian Wigner basis'^ for another 
expansion of the same G. The Wigner basis comprises the operators 

2(-i)(«^ - ^*){« - ^) = 2(-l)"^" e^^^^ , (30) 

which are obtained from the parity operator (—1)" " by unitary displacements, 
the factor of 2 normalizes them to unit trace. This gives 

G(a) = j (dz) 2(-l)('^^ - ^*)(" - ^)w^(z) , (31) 

where 

W{z) = Tr {2(-l)('^^ - ^*)('^ - ^)G(a)} 

= Tr{2(-l)'^^«G(a+z)} (32) 

is the Wigner function to G(a). A real Wigner function, [W^(z)]* = W{z), is asso- 
ciated with a Hermitian operator, [G(a)]^ = G(a). 

Since Fourier transformation relates the bases to each other, 

e^'E^ = y'(dz')e^'E^'2(-l)(-' -^'*)(«-^') , 

2(_i)(«^ - ^*)(« = j (dz') e^^^^'e-^^^^' , (33) 



''Concerning Weyl's unitary operator basis, the seminal papers by Weyl (1927) and Schwinger 
(1960) are recommended reading. ^^'^^ A recent textbook account is given in chapters 1.14-1.16 
of Ref. 21. 

"^Concerning Wigner functions, the seminal papers by Wigner (1932) and Moyal (1949) arc recom- 
mended reading, ^^'^"^ and so arc the more recent reviews by Tatarskii,^'* by Balasz and Jennings. 
and by Hillery et al.,'^'' and also the textbook expositions by Scully and Zubairy,^'^ and Sclileich.^*^ 
'Perhaps the first to note the intimate connection between the Wigner function and the parity 
operator was Royei;"^'* in the equivalent language of the Weyl quantization scheme the analogous 
observation was made a bit earlier by Grossmann.^'^ A systematic study from the viewpoint 
of operator bases is given in Ref. 31. — The appearance of the parity operator is central to 
experimental schemes for measuring Wigner functions directly. 
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the (real) Wigner function W and the characteristic function C are Fourier trans- 
forms of one another, 



(34) 



W{z)= J{dz')e--'^-'C{z'), 

C(z) = J{dz')e--^^'w{z'), 

and, therefore, the Wigner function is also a Gaussian, 

W{z) = e-^^^W^ (35) 

where 



Tr {G} = 1 is equivalent to j (dz) W{z) = 1 . 



(36) 



Note that [C(z)] * = C(-z) and [W{z)] * = W{z) imply each other. 

Perhaps the simplest verification of (33) combines the Baker-Hausdorff identity 



gztEa _ Q—za^ + z*a 



= :e^^E^-3^^^: (37) 
and the normally-ordered form of the displaced parity operator, 

e^^Ez 2(_l)«ta gztEa ^ . 2e-{a+ - zt)(a - z) . ^ ^gg) 

with the basic Fourier-Gauss integral 

(dz) e-h-'^-e^'- = ^^e^^TA-ix ^ (39^ 
vdet A 

valid for all matrices A = TA^'^T > and all columns x, whether x"^T is simply 
related to or not. 

Upon using (37) in (12) or (38) in (31) we find the normally ordered form of 
G(a), 

G(a) = \/d^ : e-^a^Qa : , (40) 

which is another Gaussian function. Fourier-Gauss integrals connect the various 
ways of writing G(a) and, accordingly, the matrices C, W, and Q must be simply 
related. Indeed, one finds 

C = EW-^E = EQ-^E- |l, 
W= EC-iE= (Q-'- 

Q = E(C + i|)-'E=(W-i + i|)-\ (41) 

where I is the 2x2 unit matrix. The symmetry property (20) of matrix C is inherited 
by matrices W and Q. 
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Since ECE, W, Q are functions of each other, these three matrices commute with 
one another. In fact, as long as we are deahng only with 2x2 matrices, Hermitian 
and with identical diagonal values, identities such as ECE = det C can be used 
to achieve a further simplification, 

C 



W = 



dctC ' 



Q 



C + i| 



(42) 



dct(C 

but this is particular to 2 x 2 matrices and does not hold for the 4x4 matrices in 
Sec. 3. 

In more explicit terms, we have 

W= (l¥{aatG} - iE)"^ , 
Q=(Tr{aatG} + i(l-E))-' 



or 



and 



W: 



Q = 



(n + If 



{n + i)2 — m*m 



n + 1 m 
m* n + 1 



4 


m 




I* n+i_ 






1 — V n 






1 - V 



(43) 
(44) 

(45) 



They obey (21 + W)(2I — Q) = 41, as they should, which one verifies easily by 
inspection. 

It is time to note that the transitions from C to W and Q are only possible if 
the respective Fourier integrals are not singular, which requires 



C > 



or, explicitly, 



n + i > V m*m 



(46) 
(47) 



Values of n and m that violate this condition will, therefore, not be considered at 
all. The determinant of C is then positive, and so are the determinants of W and 

Q, 

det C = (n + i)^ - m*m > , 



detW= (detC)-^ >0, 

detQ= I ^ >0. 

[n + 1)^ — m*m 



(48) 



2.2. Positivity criteria 

Owing to its simple Gaussian form, operator G must be unitarily equivalent to the 
basic Gaussian Gq, 



G'o = (l-ff)5"'" = (l-5) :e-(l-^)»'»: 



(49) 
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where —l<g<l ensures a finite trace. In other words, we have 

G = U^'GoU 



(50) 



with some unitary U that effects a linear transformation on a and , a squeezing 
transformation in the jargon of quantum optics. Its most general form is 



U''aU = ac''^ cosh 9 + ah''^ sinh 9 , 
U^'a^'U = ah-^'f' cosh (9 + ae"'"^ sinh 61 



or, compactly, 

U''aU = U'' 

with 



= Ua, [/tatC/ = a^U^ 



U = 



e''Pcosh6' e'"^ sinh 6* 
e-'"^ sinh 61 e"'"^ cosh 61 



(51) 



(52) 



(53) 



which is characterized by three real parameters: 0, c6, ip. In the present context, only 
the relative phase (j) — f enters, so that the initial parameters n and m determine 
g, 9, and (j) — ^p- 

Note that the 2x2 matrix U that is thus associated with the unitary operator 
U is not a unitary matrix itself. Rather it obeys 



U^EU = E or U"^ = EU^E 
to maintain the fundamental commutation relation 

aa^ — a^a = 1 or a^Ea = —1 . 

and in addition 



(54) 



(55) 



(56) 



must hold for consistency with 



atT. 



The resulting relation between the characteristic functions of Go and G amounts 



to 

C = U+CoU 

and, as a consequence of (41) in conjunction with (54), we find 

W = U^WoU 

and 



Q 



(UtQoU) ' + iE(l-UtU)E 



(57) 
(58) 
(59) 
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We remark that the statement (58) about the Wigner functions is, of course, con- 
sistent with the general observation^ that Hnear transformations on the operators 
a, are reflected by exactly the same transformation on 2, z* in W . 

Further, we note that Q 7^ U^QqU unless U commutes with E, and this is as it 
should be: transformations with EUE 7^ U turn a into a linear combination of a and 
a\ so that the meaning of normal ordering is altered. The extra term in (59) takes 
just that into account. 

Let us use the matrices C and Co of the characteristic functions to find the 
relations between the initial parameters n, m and the new parameters g,6,4> — (p. 
With C of (16) and 



in (57) we have 



and, in particular, 



Thus g is given by 



m=l li^e-'(<^ - ^) sinh(2^) , (61) 

2 1 Q 



{n+\f-m*m=\[^^^] . (62) 



^Ain^\f-m*^Y-\^ (63) 
[{n+\Y-m*mY +5 

and since the eigenvalues of G arc (1 — g)g^ with fc = 0, 1, 2, . . . wc find that G > 
requires 5 > 0, which in turn says that the argument of the square root must be at 
least \, and this is precisely the constraint (25). In other words: Condition (26) is 
both necessary and sufficient for the positivity of G. 

These considerations are instructive, but they arc not really needed if wc just 
want to find the eigenvalues of G, for which purpose knowledge of 9 and ((> — (p \s 
obsolete. A direct method proceeds from the observation that 

00 -. 

fe=o ^ 



SThc statement is actually true for rather arbitrary linear similarity transformations, not just for 
linear unitary transformations, and it applies to iimltidimensional Wigner functions. Somewhat 
surprisingly, this important transformation property is not as widely known as it should be. 
Various special cases are demonstrated in Refs. 35, 31, and 36: lineaj" unitary transformations 
for one degree of freedom; linear similarity transformations (unitary or not) for one degree of 
freedom;*^ lineaj: unitary transformations for many degrees of freedom.^® And Ref. 37 deals with 
the general case of linear similarity transformations (unitary or not) for many degrees of freedom. 
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and employs the Wigner function for calculating this trace, 

Tr{G^}= J (dz) [W{z)]^ = i Vd^ . (65) 

In conjunction with (48) this reproduces (63). 

Note that condition (65) follows from the general property of a density operator 
for which 

Tr {G^} < 1 is equivalent to j (dz) Wizf < 1 . (66) 

The discussion above tells us that condition (66), which leads to 

VdetW < 2 , (67) 

is both necessary and sufficient for the positivity of G in one dimension. As we 
shall see in Sec. 3, however, this condition is not strong enough to guarantee the 
positivity of G in two and more dimensions. 

We do not even need to know the eigenvalues of G if checking G > is all that 
we are interested in. For, the normally ordered form (40) reads more explicitly 

G = v/d^e-s'^"^' : e'^^ " '')«^« : e-^"*"^ . (68) 
This has the structure G = S^GS with some S and 

G =(!-!/) : e-(^ - : = (1 - iy)!/"^" , (69) 

or 



G = ^JdetQ : e-^^^'^^ : (70) 



with 



Q = 



1-1/ 
1-1/ 



^;(Q + EQE). (71) 



2' 

Consequently, G > is ensured by G > 0, and this just requires v >0, ov 

I - Q > . (72) 
Now, upon recalling how u is related to n and m in (45), 

_ n{n + l)-m*m 

(n + 1)^ — m*m, 

we find, once more, that (26) is the positivity criterion. 

Having found the value of <?, the unitary transformation of (50) can be identified. 
For this purpose we return to (57) and write it in the equivalent forms 

CEUt = U^CoE and UEC = ECqU , (74) 

where 



EC -11 + ^ 



1 
-1 



CoE. (75) 
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Therefore, the eigenvalues of the 2x2 matrices EC and CE are ±^(1 + — g), 
the columns of are the respective cigencolumns of CE, and the rows of U are 
the eigenrows of EC. With (61) and (63), it is a matter of inspection to veriiy these 
statements for C of (16) and U of (53). 



2.3. P-representable positive Gaussian operators 

For g = 0, the basic Gaussian Go of (49) is : e~"^" :, the projector to the oscillator 
ground state. Unitary displacements turn it into : e~("^ ~ ~ ^) which project 
onto the coherent states, the eigenbras of and eigenkets of a with respective 
eigenvalues z* and z. 

A positive Gaussian operator, G > 0, is said to be P-representable if one can 
write it as a mixture of coherent states, 

G = J{dz) :e-(«^-^*)(«-^) : Vd^e'^^^^^ (76) 

with P > 0. The limiting case of G =: e~"^" :, when'' 

Vd^e-^^''^-^<5(z), (77) 

so that P — » 00 in some sense, need not concern us too much. For the sake of 
notational simplicity, we exclude it by considering only G > 0, rather than G > 0. 
For a P-representable Gaussian, we have 

Q = (P-^ + !)-', 
w=(p-i + ii)-\ 

C = EP-iE+i|, (78) 

and 

P = E(C - i|)"'E = (W-^ - ^\y' = (Q-1 - I)-' . (79) 
So, a given G is P-rcprcscntablc if 

C-i|>0, W-i-i|>0, Q-i-l>0, (80) 

which requires 

n> \m\. (81) 

There are, therefore, positive Gaussians operators that are not P-representable, 
namely those with 



m > n> 



V"^*rn+i-|. (82) 



''See (101) below for the definition of the two-dimensional Dirac delta function. 
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All positive Gaussian operators are, however, unitarily equivalent to a P-repre- 
sentable one, because Go of (49) is P-representable if 5 > 0, as 

(1 - 5)5»^« = J (dz) : e-(«^ - - ^) : {g-^ - l)e-(«"' " (83) 

shows explicitly. If G > is P-representable, then the unitary transformation of 
(49) amounts to 

(84) 



(U^PoU) ^ - iE(l- Utu)E' 



where we see an extra term, similar to the one in (59), which reflects the injunction 
of normal ordering that is inherent in (76). 

The unitary transformation of (50) that relates the given positive G to the spe- 
cial P-representable Go of (49) and (83), is just one transformation of many, which 
are all such that UGW is a P-representable Gaussian. In terms of the respective 
characteristic functions, these transformations are such that 



C = 



n + 
m* 



n + 



with 



Ut ^CU"^ = EUECEU^E 



n>\m\, 



N- 



M 



M* -h 



(85) 



(86) 



where U is of the form (53) and obeys (54). For example, with Re(me''^ = |m| 
ensured by choosing <j) and ip accordingly, we get 

AT + i = (n i) cosh(2(9) - |m| sinh(26l) , 

M = e''f' + [|m| cosh(26l) - (n + i) sinh(26l)] , (87) 

and then (86), the condition that UGW is a P-representable Gaussian, amounts to 

1 - ,26* 



2n -h 1 - 2 m 



< e^" < 2n+l + 2 TO . 



(88) 



Owing to (47), the lower boimd is assuredly positive, and since (25) holds for a 
positive G, the upper bound is certainly larger than the lower one so that the range 
for e^^ is not empty. 

Any 6 from this range will serve the purpose of relating the given G to a P- 
representable one. In particular we note that = is permissible (of course) if G 
itself is P-representable (n > |to|). The special transformation of (50), for which 
UGW = Go and thus M = 0, obtains when e^^ is the geometric mean of the 
bounds, so that 



n-l 


4 + 


|to| 


n-l 


1 
2 


TO 



(89) 



identifies the 6 value for this distinguished squeezing transformation. 
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2.4. Transposed Gaussians 

As a preparation for a later discussion in the context of Gaussian states of two 
entangled harmonic oscillators, let us briefly discuss what happens when an operator 
transposition is done. First of all, we must note that transposing an operator is a 
representation-dependent operation. For a chosen basis of Hilbert state vectors, \k), 
the transpose of G is defined by 

{k\G^\k') = {k'\G\k) . (90) 

Clearly, if the |A;)'s are eigenkets of G, there is no difference between G and G^. 

We shall consider two procedures based on the position and momentum repre- 
sentations associated with the eigenstates of (j = 2^2 (q/^ + a) and p = 2~2 (ia'l' — ia) , 
respectively. Transposition in the ^-representation has this effect on products of a 
g-function and a p-function, 

iT 



[.Mq)f2{p)V =f2{-p)fl{q): 

and in the p-representation one gets 

[/l(9)/2(p)]^ = /2(p)/l(-g). 

it 



We find the resulting transformation of e^ 



by first noting that 



^za< — z a 



^2 2(z — z*)q-2 2i(z + z*)p 

J2-i(z - z*)q^-2-ii{z + z*)p^-^{z - z*){z + 2*) 



(91) 
(92) 

(93) 



so that 



^za' — z a 



^ Q2-h{z + z*)p^2-i{z - 2*)gg-i(z - z*){z + z*) 
^ ^2-i{z - z*)q + 2-ii{z + z*)p 



= e 



za — z a 



(94) 



in the ^-representation, whereas 



zat — z*a 



^ Q-'i-hiz + z*)p^-2-h(z - z*)q^-l{z - z*){z + z*) 
_ ^-2~^{z- z*)q-2~^\{z\z*)p 
z*a^ — za 



= e 

in the p-representation. These are compactly summarized in 

^ gTztTEa ^ giztETa 



(95) 



gatEz- 


T 


g-ztEa- 









(96) 



where the upper sign refers to the g-representation, and the lower sign to the p- 
representation. Not surprisingly, we encounter the transposition matrix T of (18). 
In G of (15) and W of (35) we thus have 

zt^iztT, 



±Tz, 



(97) 
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and 

a ^ TTa , ^ Ta^T (98) 

apply in the normaUy ordered form of (40). Since the sign is irrelevant for the 
Gaussians that wc arc concerned with, we have in both cases 

C^TCT, W^TWT, Q^TQT, P^TPT, (99) 

which arc, of course, consistent with (41) as well as (78) and (79). 

Indeed, the net effect is simply m ^ m* , and therefore has the same eigen- 
values as G, so that G and G'^ are unitarily equivalent. The situation is markedly 
different when a partial transposition is done on a Gaussian state of two harmonic 
oscillators. 



2.5. Examples 

We conclude this section on one-dimensional Gaussian operators by giving two 
explicit examples. 

2.5.1. Parity Operator 

The first example is the parity operator that is used to form the Wigner basis. 
We have previously noted that the transition from C to W and Q is nonsingular 
if C > 0. Actually, the limit C > 0, C can be included with a bit of caution, 
and C = is actually needed in the construction of the Wigner representation 
of the Gaussian operators (31). For this choice of C, we obtain: Q = ^1, and 
accordingly the corresponding Wigner function of such an operator is given by a 
singular expression that can be normalized: 

W(z) = 5{z) with j (dz) W{z) = 1 , (100) 

where 

^(z) = 2n6iq')dip') , J (dz) <5(z)/(z) = / (z = 

for both parameterizations (8) and (14). Using this relation, we obtain from (31) 
that the resulting Gaussian operator 

G = 2 :e-2«^": = 2 (-1)"^" (102) 

is twice the well known parity operator. This operator is non-positive and normal- 
ized to unit trace, Tr {G} = 1. Note that the Wigner function (100) is singular, 
normalized and positive, while the corresponding Gaussian operator is normalized 
but not positive. Upon unitarily shifting the parity operator by complex numbers 
z, we reproduce all elements of the Wigner basis (30). 



) (101) 
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2.5.2. Pure Gaussian State 

The second example is a Gaussian operator describing a pure quantum state. It 
follows from (65) that a Gaussian state is pure if 



detC=(detW)-i = (n+i)2-m*m=^ (103) 
holds, which amounts to 



|m| = V"(n + 1). (104) 
In this case the general formula (68) reduces to a projector 



G = = Vd^e-2^"'^ |0)(0| e-2/^*« , (105) 



where 



and 



detQ = - = 1-11*11, (106) 



n+ 1 ' n+ 1 



|0)(0| =:e-"^" : (107) 

projects on the 0th Fock state, the ground state of the harmonic oscillator: a|0) =0. 

In (104) we recognize the border case of (25), and thus of (24). This tells us 
that one can check the purity of a Gaussian sate of a one-dimensional oscillator by 
measuring the expectation values of a^a and a?, and then verifying that the equal 
sign holds in (24). 

From (105)-(107) we see that the ket vector of the pure Gaussian state has the 
form 

= (l-///i)ic-i^"^'|0). (108) 

Now assuming, for simplicity, that /x is real and positive, /x = > 0, we recognize 
in this expression a squeezed state of a one-dimensional harmonic oscillator, whose 
position wave function is 

V'M(9) = (9|M) = (^)'e-W (109) 

with 



K = = , — . (110) 

1-M y/nTT-^ 

For n = 0, we have ii = 0, k= 1; the wave function is then that of the ground state 
of the one-dimcnsional harmonic oscillator, as it should be. 

The general pure squeezed state is never P-representable because the relation 
(81) is never satisfied when (104) holds. But as we have seen above, it can be 
unitarily transformed into the projector |0)(o| which is a P-representable Gaussian 
operator. 
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3. Gaussian States of a Two-Dimensional Oscillator 
3.1. Pararaeterizations 

Things look much the same except that the number of variables doubles. Thus now 
we have a\, ai and zl, z\ for the first degree of freedom, as well as ai and z^-, zi 
for the second, so that 



ai 



at = 



a|, ai, a|, 



z\ 
z*i 

Zl 

Vzl 



= [zl,Zi,Z2,Z2\ 



(111) 



are 4-component columns and rows, and C, W, Q, P are 4 x 4-matrices. The number 
of integration variables doubles as well, of course, and 



J{dz)...= J{dz,) j{dZ2)... 



(112) 



specifies how phase space integrals are to be understood. Relations (41), (78), and 
(79) remain valid with 



10 
0-100 
10 
-1 



(113) 



replacing the 2x2 version of E in (11). 

Linear unitary transformations of the form (50) now involve a 4 x 4 matrix U 
for the transformation of a and a^' as in (52), which continues to obey (54) and (56) 
with 



10 
10 
000 1 
10 



(114) 



now. With these replacements, the 4x4 matrices C, W, Q, and P transform ac- 
cording to (57), (58), (59), and (84), respectively. 

The explicit parameterization of the matrix appearing in the Gaussian charac- 
teristic function is 



(115) 



"ni + i 


mi 


rris 


rUc 


m\ 




m* 


m* 


m* 


rric 




1712 




rris 


m2 


n2 + 



and C > is taken for granted again. Here, too, cxpaiiding 



e 2^ 



Tr {e^^^^c} = Tr {f 



(116) 
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in powers of z reveals the physical significance of the numerical parameters in C, 
viz. 



m 


= Tr • 


a[a^G 


mi = 


-Tr {alG} , 


n2 


= Tr ■ 


a\a2G 


m2 = 


-Tr {alG} , 


rus 


= Tr . 




TOc = 


— Tr |aj^a2G' 



(117) 

and positive two-dimensional Gaussians must obey (26) for analogous reasons as in 
the one-dimensional case of Sec. 2. 

Partial transposition, in the 21,^1 sector only, turns G into G'^^ . For C, W, Q, 
and P the transition is as given in (99) with T replaced by 



Ti = 



10 
10 
010 
1 



(118) 



We note that 









mi - 


-> —mi 


c- 


^ ECE 


amounts to < 


m2 - 


-m2 








. mc - 


-rUc 



(119) 



with ni,n2,ms unaffected, and 









' mi - 


mi 


c- 


TiCTi 


amounts to < 


nis - 










, rric - 


TO* 



(120) 



with ni,n2,m2 unaffected; analogous statements apply to Q ^ EQE, Q TiQTi, 
et cetera. 



3.2. Positivity criteria 

A first positivity check is the one that exploits that G must be unitarily equivalent 
to a (product of) basic Gaussian(s), 

G = C/t (1 - 5i)5^^"i (1 - 52)52^"' U (121) 

with some appropriate unitary U. The positivity of G can be tested by checking 
whether gi and g2 are positive, for which it suffices to see if the smaller one of the 
two is positive.' 



'In marked contrast to the one-dimensional case, condition (66) alone does not guarantee the 
positivity of a two-dimensional Gaussian operator. As a counter example put 91 = 4, 92 = —5 in 
(121), so that the resulting Gaussian has negative eigenvalues while Tr {G} = Tr {G^} = 1- We 
are reminded here that the positivity of the statistical operator for a two-party system requires 
in particular the positivity of the reduced density operators, here characterized respectively by gi 
and 92 ■ In two dimensions, we also need to consider Tr { } when verifying the positivity of G. 
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Now, information about gi and 52 is available in the traces 
TrjG^} 



^21 _ 1 -51 1 -52 



1 + 51 1 + 52 ' 

Tr{Gn = t^i\^. (122) 



The trace of G can be expressed in terms of the C matrix, 



Tr {GH = = , (123) 

^ ^ V^det(2C) 4VditC 



and upon introducing 

i2 



G ^ = C/t (1 - 5?)5^'"^ (1 - 52^)52"^"^ : 



(124) 



with matrix C in its characteristic function (see Sec. 3.3 below), we have 



(Tr{G2})2 I J 4Vd;^- 
So, we extract the necessary information about gi and 52 out of j 

4Vd^=i±^i±^ and 4v^= ;+i;+i . 

1 - 5i 1 - 52 1 - 5f 1 - 5| 



(126) 



A straightforward, yet somewhat tedious, calculation then establishes that G > 
holds if both inequalities in 



+ 3 dot C > 4VdetC VdetC < i + 2 det C (127) 

are obeyed, and only then. In deriving these two inequalities, it is taken for granted 
that 

-1< .91,52 <1, (128) 

which is not an actual limitation because G^ = G and Tr {G^} < 00 restrict the 
values of 51 and 52 to this range. 

The right inequality in (127) amounts to 

5152 > (129) 

and the left to 

(51 + 52)(1 + 5i)(l + 52) > (51 - 52)' . (130) 

The first requires that both 5's are positive or both negative, but the second excludes 
the possibilities 51 < 0,52 < and 51 < 0,52 < 0. So, taken together, the two 
inequalities of (127) imply 51,52 > 0, that is: G > 0, indeed. 



■•Once gi and 92 are available, a variety of unitarily invariant quantities can be computed, such 
as the von Neumann entropy of the mixed state represented by G.^^ 
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Alternatively, for a second positivity check we can apply the reasoning of (68)- 
(73) to 



Q 



\-Vi fil Us 

/i* Us M2 1-1^2. 



(131) 



where now 



Q 



1 — Pl jJLs 

/u* 

/i* 1-1^2 

lis 1 — 



2(Q + EQE) 



(132) 



shows up in the 4x4 version of (70). Since Q commutes with E, we only need to 
consider unitary transformations with UE = EU and thus U^'U = I when diagonal- 
izing G, so that the complication of the extra term in (59) is of no concern here. 

Therefore, the criterion (72) applies to the two-dimensional case as well, and we arc 
asked to check if the eigenvalues of Q exceed unity or not. These eigenvalues are 



1 - + ^2) ± \l[\{yi-l^2)]^ +tl>s 



so that G > if 



or 



and only then. 



J^l + J^2> V {Vl - V2y + 4/U*/Us 

yi + V2>^ and viV2 > nlUg 



(133) 
(134) 
(135) 



3.3. Squared Gaussians 

The positivity criterion (127) needs knowledge of C, the 4 x 4-matrix appearing in 



the characteristic function of G. It is found by squaring 

G = y^(dz)e"^^^-5^^^^ 
recalling the Baker-Hausdorff relation et cetera, with the outcome 



C=iC + iEC-iE=iC+iW. 



As a consequence, then 

detC= j;^det(C + iEC-^E) , 
and (127) reads more explicitly 



j^+3detC> A/(detC)det(C-f- iEC-^E) < i-h2detC. 



(136) 

(137) 
(138) 

(139) 
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3.4. Separable and non-separable Gaussians 



As a rule, the statistical operator p(^'^^)(a) 



,(1&2) 



(31,82) of a two-dimensional 



oscillator, be it of Gaussian shape or not, will be different from the product of the 
reduced statistical operators p^^^(ai), p^^^(a2) that one obtains by partial tracing, 



with 



p(^^^)(a)^pW(a0p(^)(a2) 



p(i)(ai) = 1^2 {p(^^^)(a)} , p(2)(a2) = TVi {p^'^^){^)] 



(140) 



(141) 



But, as we mentioned in the Introduction, it could happen quite easily that 

is the convex sum of such products, in which case it represents a separable state, 

p^^^^) is separable if p^^^^^ = ^ Wkp^j^^ p^^'^ with Wfe > . (142) 

k 

Whereas it can be rather difficult to decide if a given pf^^^) jg separable, matters 

are remarkably simple for Gaussian states. 

Positive Gaussians that are separable must have a positive partial transpose 

this is Pcrcs's necessary criterion. In the case of Gaussian states, it is also 
sufficicnt.14.15.40 

Concerning the positivity of G^^ , please note that dct(TiCTi) = det C and that 

1\a\ = iJi (C + iTiETiC-iJiETi) Ti 
= iTi (C + iET^C-iET^)Ti 

so that a non-negative partial transpose, G'^^ > 0, is available only if 



^-h3detC> ^(detC)det(C-hiETiC-iETi) < i-h2detC 



(143) 
(144) 



hold with 



TiETi 



-10 
10 
10 
-1 



(145) 



Probably simpler is the positivity criterion of (131)-(135). When applied to G'^^ 
it asks whether the eigenvalues of 



-(TiQTi + ETiQTiE) 



n* 
1 — i/1 pc 

Pc 1 — 1^2 

Pi I-U2 



(146) 



exceed unity. Therefore, if G > so that the restrictions of (135) are obeyed, G'^^ 
is positive if 



Z^l!^2 > P*cPc ■ 



(147) 
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Gaussians that are P-representable, 



(148) 



(with P > 0), are separable by construction. Local unitary transformations, Uio 
U1U2 with 



Ui 



and 



e'?"! cosh fill e'^'isinhfi'i 

c-'^isinh^i e-''^icoshei 

10 

1 



10 

10 
6*^2 cosh 6I2 e'<^2sinh6l2 
e-''^2sinh6i2 e-'-^a cosh 



(149) 



(150) 



turn such a P-representable Gaussian into other separable Gaussians, which may or 
may not be P-representable themselves. All separable Gaussians can be constructed 
in this way, and as a consequence the Peres criterion is not only necessary, it is 
indeed sufficient. 

Although this is reasonably obvious, the explicit construction of the unitary 
transformation that does the job is technically more involved than the analogous 
problem in one dimension that we treated in Sec. 2.3. We shall, therefore be content 
with the following general remarks, and present an explicit example in Sec. 3.5.4 
below. 

The separability of a Gaussian state can be investigated by studying the ex- 
istence of a local transformation U\oc = U1U2 that maps the general Gaussian 
characteristic function (115) into a P-representable Gaussian characteristic func- 
tion 



t^ioc : C 

which has to satisfy the condition 



Cp- i| >0. 



(151) 



(152) 



This says that the four eigenvalues of Cp — ^1 are strictly positive. As shown in 
Refs. 14 and 15, this is indeed equivalent to the separability condition (5). 



3.5. Examples 

We conclude this section on two-dimensional Gaussian operators by presenting a 
couple of examples, mostly generalizations of the historical EPR Gaussian state. 
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3.5.1. Pure Gaussian state 

We begin our examples with the simplest case of a pure Gaussian state of a two- 
dimensional harmonic oscillator. For pure states, we have Tr{G^} = 1, and so 
(123) implies that the Gaussian characteristic function must have a C matrix (115) 
with 



detC : 



16 



For convenience, we rewrite (115) in the block form 



C = 



Cll Ci2 

C21 C22 



(153) 



(154) 



For such pure states, the criterion for separability is simple, inasmuch as it is 
enough to study the properties of Cn or C22, which correspond to the reduced 

one-dimensional parts of the two-dimensional harmonic oscillator. If both Cn or 
C22 describe pure states, we say that the joint state is not entangled. This amounts 
to 



det Cll 



1 



or det C22 



(155) 



For Gaussian pure states this separability condition is much simpler than the general 
condition of Sec. 3.4 for mixed Gaussian states. 

If we restrict the parameters that characterize a pure Gaussian state to real 
numbers only, and continue to ignore linear shifts, such a state is described by a 
wave function of the form 



*(q) 



dctD\ ^ 



(156) 



in the position representation (q = [gi,(j'2]^)- Here we meet the 2x2 matrix 



D = DT 
and D > 0, that is 



a 7 
7/3 



with det D = a/3 - 7^ 



a + /3> y(^^^^pT4^, 



(157) 



(158) 



is needed for the normalizability of the Gaussian wave function (156). 

Then, the C matrix of the resulting Gaussian characteristic function (154) has 
the blocks 



-11 



/3 



a p a 

4 ^ 4detD 4 4detD 

a _ a , (3 
4 



a 

4detD 4'^4detD 
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-22 



/3 a 
4 4detD 
(3 _ a 
4 



4 



a 



4detD 



Cl2 = Coi 



1 - 

1 + 



4detD 
1 



4 4detD 



1 



dctD 
1 

detD 



1- 



detD 
1 

detD 



(159) 



It is easy to check that a C matrix with these blocks satisfies the purity condition 
(153). 

The separabihty condition det Cn = \ amounts to det D = a(3, so that the state 
is separable (not entangled) if and only if 7 = 0, as it should be expected for pure 
states. For 7 7^ 0, all pure states represented by the Gaussian wave function (156) 
are entangled. 

This wave function defines a Gaussian operator that is a projector, 



where 



with 



G= |^')(«'| , 
I*) = (detQ)3e-5^'i"l'e-|/^24'e-Mc44 |o^o> 
det D + a - /3 - 1 



(160) 

(161) 



Ml = 



M2 = 



det D + a + /3 + 1 

det D - a + /3 - 1 
det D + a + /3 + 1 

27 



det D + a + /3 
being the non-zero parameters in (131), and 

16 det D 



1 



detQ = 



(162) 



(163) 



(det D + a + li + lY 

follows. Therefore, such a Gaussian wave function corresponds to two squeezed one- 
dimensional harmonic oscillators, correlated with a strength that is characterized 
by the parameter 7. 

As a further simplification, and also as an important physical application related 

to the EPR wave function (2), let us investigate a correlated state that is not 
squeezed. We obtain such a state by choosing the parameters in accordance with 

a = (3 and det D = 1 , so that -f'^^a^-l. (164) 

Such a Gaussian state is characterized by a single real parameter, for which n in 

a = (3=l + 2n and I7I = 2^Jn{n + l) (165) 
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is a convenient choice. Upon opting for 7 = — I7I, the wave function (156) is 

*(q) =^-^ cxp(^-{n+^)iqf + ql) + 2^nin + l)qiq2) . (166) 
For large n it reduces to 

*(q) ~e-^(^i-''2)2 (167) 

and becomes 

*(q) ~ S{qi - q2) (168) 

in the hmit n 00. We recognize here the famous historical wave function (2) of 
two correlated particles introduced by the EPR trio. This wave function cannot be 
normalized, and thus represents an overidealized situation. By contrast, the wave 
function (166) is a smoothed, normalized version of (2) and refers to a real physical 
situation. For illustration, we show this smoothed wave function in Fig. 1 for two 
different values of n. 




"1 % 

Fig. 1. Contour plot of the smoothed EPR wave function (166) for n = (not entangled, on the 
left) and n = I (entangled, on the right). 



3.5.2. Bell states 

For entangled spin-i states, a special role is played by the entangled basis of the four 
Bell states. For the two-dimensional harmonic oscillator it is possible to introduce 
a continuous- variable generalization of Bell's states. Such a complete set of states 
consists of particular Gaussian operators. 

Following the example (100) of the singular Wigner function for a one- 
dimensional harmonic oscillator, we write for the two-dimensional case a singu- 
lar Wigner function in the form of a product of two two-dimensional Dirac delta 
functions. 



W{z) = 6{zi - zo)(5(z2 + zo) with z = 



and zq = 



(169) 
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where the two-dimensional (5-function 6{z) is as defined in (101) and zq is an arbi- 
trary complex number. The corresponding Gaussian operator, 

G{zq) = • Q^'*!"! ~ "2"2 + "l"2 -I- a}a2 + ^0(12 - al) -I- ^5('*2 ~ "l) • 

=|zo)(zo|, (170) 
projects on the "continuous Bell state" specified by the ket 

|zo) = Q-WM'^ + -4A-zo<A + 'AA |o,o) 

= eil^ol^ + ("1 - + ^0) 0) , (171) 

which obtains from |zo = [0,0]""") by a simultaneous unitary shift of ai by zq and 
32 by — z^, respectively. The completeness relation 



'(dzo)|zo)(zo| =1 

is easily demonstrated with the aid of the Fourier-Gauss formula (39). 



(172) 



3.5.3. Mixed EPR states 

In this example the matric C of the Gaussian characteristic function is of the form 



n + 













rric n -J- 



n + i 



rric 









n-h 5 



(173) 



so that n (= ni = 712) and rric are the only non-zero parameters in (115) and, in 
addition, we take rric to be real. The corresponding Q matrix of (131) is 



Q = 



{n + iy 



mf. 



n -|- 1 rric 
n-\-l ruc 
rric n-\-l 
rri c n -M 



and the positivity criterion (135) amounts to 

n{n -f 1) > . 
This is less stringent than the condition for separability, 

n > \mc\ , 



(174) 



(175) 



(176) 



that follows from (147). These matters arc illustrated in Fig. 2. 

Note that the Gaussian operator specified by (173) represents a mixed state in 
general. Only if 

1 

16 



detC = [(n+ 5)^ 



212 



(177) 
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2i 




Fig. 2. Criteria for positivity and separability for the mixed EPR state specified by (174). In this 
plot of n versus rric, the states associated with points on and above the dashed line are separable, 
whereas points on and above the solid line mark positive Gaussians, with pure states obtaining 
for n, rric values on the solid line. 



holds, the state is a pure state, which is the border case of the positivity criterion 
(175). With the identification 7 = 2toc and n = n, we then recover the EPR pure 
state of Sec. 3.5.1, as we should. 

We further note that the remarkably simple separability condition (176) can be 
expressed in terms of Bohr's variables (4) as 

{Vl) + {Ql)>l and {Vl) + {Ql)>l. (178) 

These relations have important physical implications when the separability of the 
EPR state is put to an experimental test. It is then enough to measure the variances 
of the Bohr variables and check if the conditions (178) are satisfied. One should re- 
member, however, that (178) applies only to the very special state characterized by 
the C matrix of (173). As we will see in the next examples, Gaussian states different 
from this generalization of the EPR state will lead to a separability condition that 
will be different in form. 



3.5.4. Anti-EPR states 

We generalize the previous example by adding anti-EPR terms to the matrix C of 
the Gaussian characteristic function in (173), so that now C is of the form 



C = 



n 
rrir 







rric 

TO., 



rric 
n 



rric 
rris 




(179) 



and is characterized by the three real parameters n, TOc, and m^. According to 
(120), the partially transposed Gaussian has TOc and nis interchanged, and because 
of this property we call the to^ terms the anti-EPR contribution. The Q matrix 
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(131) of such a squeezed Gaussian operator has these non-zero elements: 

(n -I- l)[(n -I- 1)^ — — m^l msmJn -\- 1) 

vi = V2 = l- — , Ml = M2 = ^ , 

a a 

^,^-m,^ , Mc = mc^ , 180 

a a 

where d = [(n + 1)^ — (m^ -I- TOc)^][(n -t- 1)^ — {m^ — As a consequence, the 

positivity condition (135) reads here 

n{n + 1) - 2m, (n + i) + > , (181) 

and the separabiUty criterion (147) amounts to 

n{n+l)-2mc{n+\)+ml-ml>0. (182) 

For mg = 0, we recognize the particular case of the mixed EPR state discussed in 
Sec. 3.5.3. 

In Fig. 3, we depict the regions for the parameters n and mc for which the 
anti-EPR Gaussian states are non-separable, for two ratios of to mc- These 
regions are bounded by the curves defined by the positivity condition (181) and the 
separability criterion (182). Note that the right plot corresponds to to, = mc, so 
that there the region of interest has no area at all, which means that such a state 
is never non-separable. 




We conclude this example by the explicit construction of the local transfor- 
mation that maps the Gaussian state associated with (179) onto a Gaussian P- 
representable state, that is 



[/loc(e) :C->Cp(0) 



M 
Ms 

Mr 



M 
M + 
Mc 
Ms 



Ms 

Mr 



N ^ 



Mc 
Ms 
M 
M^ N + 



(183) 
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in terms of the C matrix of (179). Because we have selected real parameters in C, 
the local transformation is labeled by a single parameter 9. The transformed matrix 
Cp belongs to a P-representable Gaussian operator if the conditions 

N±M>\Mc±Ms\ (184) 

are obeyed. Here we have 

iV±M= (n+i)eT2e_ 1^ 
\M,±M,\=c^^^\mc±ms\. (185) 
The lower bounds in (184) are achieved if 

2M = \Mc + Ms\-\M,-Ms\, 

2N=\M, + M,\ + \M,-M,\. (186) 

The first equality implies 

(n + |)(e-2^ - e^^) = |mc + ms\e~^^ - \mc - ms\e^^ (187) 
or, more explicitly, 

^4. ^ n+l-\rn. + mA ^^^^^ 
n + 5 - |toc - nisi 

for the transformation parameter 9. We see that for the two special cases of the 
EPR state (mg = 0) and the anti-EPR case {rric = 0), we have ^ = 0, and no local 
transformation of the Gaussian characteristic function (179) is required. In these 
two cases the P-representability can be studied directly by using (179). 

More generally, the conditions for P-representability that follow from (184) and 
(188) are 

(n + 5) - ^e"^^ > |mc + ms\ , 
(n+ i) - ie^^ > |mc-ms| . (189) 

They are easily seen to be equivalent to the separability criterion (182). 



3.5.5. Squeezed EPR states 

As the last example of this tutorial we consider squeezed EPR states. Their C 
matrix (115) is specified by the three real and positive parameters n, rric, m and is 
of the form 



m 


rrir. 



n + 
rUc 






n + 



rric 

m 

m n + -, 



(190) 
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Note that this matrix differs from the one in (173) by the additional squeezed terms 
associated with m (= mi = m^)- The Q matrix (131) of such a squeezed Gaussian 
operator has the non-zero elements 

(n-h l)[(n + 1)2 - - m^l (n + if - + ml 
vi^i'2 = l ^ , Ail = = m — , 

\2 ^2 I ^2 



2mmc(n + 1) (n+ 1) — mf^ + m 

■■ — irtr 



D ' ^ D 



(191) 



with D = [(n + 1)2 - (to + mcf][{n + 1)^ - (m - rricf]. 
Here the positivity condition (135) requires 

n{n + 1) - {rric + mf > (192) 

and, according to the separability criterion(147), the squeezed EPR state is sepa- 
rable if 

n(n + 1) - 2mc{n + ^) + ml - m'^ > (193) 

holds. For to = 0, we recognize the particular case of the mixed EPR state discussed 
in Sec. 3.5.3. 

In Fig. 4, we depict the regions for the parameters n and TOc for which the anti- 
EPR Gaussian states are non-separable, for two ratios of mtomc- These regions are 
bounded by the curves defined by the positivity condition (192) and the separability 
criterion (193). 




I 2 I 2 



Fig. 4. The positivity and separability criteria (192) and (193) bound the region of non-separable 
squeezed EPR states (solid and dashed curves, respectively). The left plot corresponds to m = 
irric, the right plot to m = rUc- 



4. Concluding remarks 

The main goal of this tutorial is to provide, for a reader trained in Quantum Optics 
or in Quantum Information Theory, a general and detailed description of various 
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mathematical tools used in the description of one-party and two-party Gaussian 
states. For such states we presented a general parametrization of Gaussian states 
and discussed the positivity criteria and the so-called P-representability of positive 
Gaussian statistical operators. With the help of these techniques wc showed that it 
is possible to address the problem of quantum separability of two-party Gaussian 
states using different-in-form, but mathematically equivalent, methods based on the 
Hciscnbcrg uncertainty relations, partial transposition, or the P-rcprcscntability. 
The tutorial includes a number of simple, instructive examples that illustrate the 
techniques and methods. 

The one-dimensional Gaussian operators discussed in Sec. 2 have been widely 
used in various application of squeezed states of light and in the theory of op- 
tical Gaussian beams involving linear optical elements. For such one-dimensional 
problems an important issue dealt with in the literature is the degree of nonclassi- 
cality of one-mode Gaussian states of the quantum electromagnetic field. A general 
discussion of this problem has been given in Ref. 41. In a different approach the 
nonclassicality has been investigated with the aid of the Uhlmann fidelity of two 
single-mode Gaussian states. 

Owing to their simplicity, Gaussian states have found many applications in 
Quantum Information Theory. Examples of recent achievements that involve Gaus- 
sian states include important problems such as entanglement key distribution,*^ 
quantum error correction,^* entanglement distillation,*^ and quantum cloning.*^ 
The separability and distillability of two-party Gaussian states are reviewed in 
Ref. 47. The problem of quantum separability of three-party Gaussian states and 
the entanglement criteria for all bipartite Gaussian states have been addressed as 
well.*^ A number of issues related to the bound entangled Gaussian states,*^ and 
the optical generation of entangled Gaussian states''*' have also been addressed in 
the published literature. 

Our tutorial is confined to Gaussian states. The reader may wish to consult 
Ref. 51, a "primer" on general composite quantum systems involving entangled 
states in discrete, finite-dimensional Hilbert spaces. 
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